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In plants, osmotically driven flows are believed to be responsible for translocation of sugar 
in the pipe-like phloem cell network, spanning the entire length of the plant. In this paper, 
we present an experimental and theoretical study of transient osmotically driven flows 
through pipes with semipermeable walls . We extend the experimental work of Eschrich, 
Evert and Young ( Eschrich et al. 1972f ) by providing a more accurate version of their 



experiment allowing for better comparison with theory. In the experiments we measure 
the dynamics and structure of a "sugar front" , i.e. the transport and decay of a sudden 
loading of sugar in a pipe which is closed in both ends. We include measurements of 
pressure inside the membrane tube allowing us to compare the experiments directly with 
theory and, in particular, to confirm quantitatively the exponential decay of the front in 
a closed tube. In a novel setup we are able to measure the entire concentration profile 
as the sugar front moves. In contrast to previous studies we find very good agreement 
between experiment and theory. 

In the limit of low axial resistance (valid in our experiments as well as in many cases 
in plants) we show that the equations can be solved exactly by the method of character- 
istics yielding, in general, an implicit solution. Further we show that under more general 
conditions the equations of motion can be rewritten as a single integro-differential equa- 
tion, which can be readily solved numerically. The applicability of our results to plants 
is discussed and it is shown that it is probable that the pressure-flow hypothesis can 
account for short distance transport of sugar in plants. 



1. Introduction 



The translocation of sugar in plants, which takes place in the phloem sieve tubes, is 
not well underst ood on the quantitative level. The current belief, called the pressure 
flow hypothesis dNobell j 1999), is based on the pioneering work of Ernst Munch in the 
1920'ies (see eg. lMunchlll930l) . It states, that the motion in the phloem is purely passive, 
due to the osmotic pressures that build up relatively to the neighboring xylem as a 
response of loading and unloading of sugar in different parts of the plant, as shown 
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Figure 1. In plants, two separate pipe-like systems are responsible for the transport of water and 
sugar. The xylem conducts water from the roots to the shoot while the phloem conducts sugar 
and other nutrients from places of production to places of growth and storage. The mechanism 
believed to be responsible for su gar-transloca tion in the phloem, called the Munch mechanism 
or the pressure- flow hypothesis (|Nobell[i999t ). states that following: As sugar is produced via 
photosynthesis in sources it is actively loaded into the tubular phloem cells. As it enters the 
phloem, the chemical potential of the water inside is lowered compared to the surrounding tissue, 
thereby creating a net flux of water into the phloem cells. This in-flux of water in turn creates a 
bulk flow of sugar and water towards the sugar sink shown on the right, where active unloading 
takes place. As the sugar is removed, the chemical potential of the water inside the phloem is 
raised resulting in a flow of water out of the sieve element. 



in figure [TJ This mechanism is much more effective than diffusion, since the osmotic 
pressure differences caused by different sugar concentrations in the phloem create a bulk 
flow directed from large concentrations to small concentrations, in accordance with the 
basic need of the plant. Such fl ows are often called Osmotically Driven Pressure Flows 
dThompson fc Holbrook 2003a ). or Osmotically Driven Volume Flows ( Eschrich et i 



19721) . 



It is, however, not clear how well this mechanism is able to account for the sugar 
translocation in plants on the quantitative level . Since the sieve tube elements that make 
up the phloem are living cells, the picture can indeed be much more complicated. For a 
large tree it would thus seem improbable that sugar transport e.g. from leaf to root by 
this mechanism would be sufficiently efficient, and in this case active transport processes 
might play an important role. On the other hand, transport over short distances, e.g. 
locally in leaves or from a leaf to a nearby shoot might be more convincingly described 
by the pressure-flow hypothesis. In any case, we need a better understanding of such 
flows in order to decide, whether they compare sensibly with translocation in plants, 
and this is the aim of the present paper. In particular, we have chosen to concentrate 
on transient flows caused by a sudden loading of sugar. First of all this gives us the 
possibility of observing dynamical behavior which allows us to compare quantitatively 
with theory and second, we can observe such dynamics using simple boundary conditions 
(e. g. closed ends), which are easily realized in experiments. 

1.1. Previous experimental work 



To study the osmotically driven flows, Eschrich, Evert and Young (jEschrich et aLlll972h 
conducted simple model experiments. Their setup consisted of a semipermeable mem- 
brane tube submerged in a water reservoir, modeling a phloem sieve element and the 
surrounding water-filled tissue. At one end of the tube a solution of sugar, water and 
dye was introduced to mimic the sudden loading of sugar into a phloem sieve element. 
The motion of this "sugar front" was monitored for different configurations of the tube 
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Figure 2. Experimental data from (Eschrich et al. 197^). Data points marked with an a repre- 
sents results from closed tube experiments and points marked with a b represents results from 
semi-closed experiments. t^ p was taken from the original paper, figures 8 and 9. To calculate 
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by observing the motion of the dyed front. They conducted experiments with the tube 
closed at both ends (closed) and open at one end and closed at the other (semi-closed). 

In the case of the closed tube, they found that the sugar front velocity decayed expo- 
nentially as it approached the far end of the tube. Also, they found the initial velocity of 
the sugar front to be proportional to concentration of the sugar solution. Through simple 
conservation arguments, which we shall go through briefly below, they showed that for 
a flow driven according to the pressure-flow hypothesis, the velocity of the sugar front is 
given by 

U/ = £ exp (~£) where to = ^Tr (L1) 

where t is time, L is the length and r is the radius of the tube, L p is the permeability of 
the membrane and II is the osmotic pressure of the sugar solution, equal to RTc (where 
R is the gas constant, T the absolute temperat ure and c the concentrat ion in moles pr. 
volume) for dilute solutions of ideal molecules ( Landau fc LifshhQ 1980h . If one applies 
this result to the flow inside a single sieve element (L = 1mm), one gets a characteristic 
velocity of ~ 7mh _1 , almost an order of magnitude larger than that observed in plants 
(see table 13]). The fact that equation [TTT] predicts a larger velocity is not surprising. One 
should keep in mind that the sieve cells consist of sieve elements separated by sieve plates 
and the resistance in strongly concentrated in the latter. Thus our experiments would 
model only the transport in a single sieve element without significantly changing the 
resistance. Secondly, the measurements in the table are supposed to be representative of 
a steady state situation and not the movement of a sugar front. 

In their experiments Eschrich et al. found good qualitative agreement between their 
results and the prediction made by equation ll.il However, as can be seen in figure the 
quantitative agreement between experiments and theory was extremely poor, the theory 



4 



K. H. Jensen, E. Rio, R. Hansen, C. Clanet and T. Bohr 



generally predicting decay times at least an order of magnitude smaller than observed. 
This discrepancy, however, is not surprising given the fact that Eschrich et al. did neither 
determine the fundamental properties of the membrane (permeability, elastic modulus), 
nor of the sugar (osmotic strength) independently. Also they did not take into account 
the unstirre d concentration layers which may oc cur near the membrane walls, effectively 
lowering IT (|Pedlevlll980l . Il98ll [l983l ; lAldislll988T) . The disagreement between theory and 
experiment, and the fact that the tracking of the sugar front was done indirectly, leave 
the details of the observed process unclear. In particular, the experiments of Eschrich 
et al. did not include a continuous monitoring of the pressure. They state (p. 295) that 
the turgor pressure rapidly builds up to a constant value, as predicted theoretically. They 
observe, however, that the tubes start leaking after around 100-150 min. (typical running 
time of the experiments), and this seems to indicate that the pressure continues to grow. 
Another source of error could come from the fact that Eschrich et al. used a dye for 
tracking the sugar front instead of directly monitoring the sugar concentration. Finally, 
it is not clear whether the membranes arc sufficiently impermeable to the sugar (sucrose) 
used. 

To make progress on these issues we have refined the experiments done by Eschrich 
et al. to test the pressure-flow hypothesis more accurately. We have measured the per- 
meability of the membrane tube, and the osmotic strength of the sugar solutions inde- 
pendently and we continuously monitor the pressure during the experiment. Also, we 
have measured the elastic modulus of the membrane tube, to asses the importance of 
elastic effects in our system. For more detailed comparison with theory it is important 
to be able to assess the entire concentration profile, and to do this, we have introduced a 
new type of experiment, where the sugar concentration is determined directly by optical 
refraction. 



2. Experimental setup 

In our experiments, we used two setups. The first setup, from now on called setup I, is 
based on the design of Eschrich et al. , but includes continuous pressure measurements. 
Further, to avoid leaking of sugar across the membrane, we use a sugar {dextran) with a 
much larger molecular weight than sucrose. The second setup, Setup II, was built in the 
spae of a prism to be able to detect the evolution of the entire sugar profile (again using 
dextran) rather than just the front position offered by previous experiments. The two 
setups are discussed in detail below. Details of the sugar and semipermeable membranes 
can be found in appendix [Al 

2.1. Setup I 

Setup I, shown on the left in figure 02 consisted of a 30 cm long, 30 mm wide glass tube 
in which a semipermeable membrane tube of equal length and a diameter of 10 mm was 
inserted. At one end, the membrane tube was fitted over a glass stopcock equipped with 
a rubber stoppper. At the other end, the membrane tube was fitted over a brass cylinder 
equipped with holder to accommodate a pressure transducer for measuring the pressure 
inside the membrane tube. 

After filling the 30 mm wide glass tube with water, water was pressed into the semiper- 
meable tube with the syringe. Care was taken that no air bubbles were stuck inside the 
tube. For introducing the sugar solution into the tube, a syringe was filled with the 
solution and then attached to the lower end of the stopcock which was kept closed. Af- 
ter fitting the syringe, the stopcock was opened and the syringe piston was very slowly 
pressed in, until a suitable part of the tube had been filled with the solution. Care was 



Osmotically driven pipe flows and their relation to sugar transport in plants 5 




i 



Figure 3. Left: Setup I used to observe the movement of a sugar-dye solution (ss) inside a 
semipermeable membrane tube (spm). L length of membrane tube; / initial sugar front height; 
ds disposable syringe; gt glass tube; rs rubber stopper; sc stopcock; wr water reservoir; be brass 
cylinder; pt pressure transducer. 
Right: Setup II. See details in the text. 

taken to avoid any mixing between the sugar solution and the water already present in 
the semipermeable tube. To track the movement of the sugar solution it was mixed with 
a red dye and data was recorded by taking pictures of the membrane tube at intervals 
of 15 minutes using a digital camera. Details of how the motion of the sugar front was 
derived from the images is discussed in appendix lB.il 

2.2. Setup II 

Setup II. shown on the right in figure [3] consisted of a hollow isosceles glass prism and 
a Plcxiglas cuboid in osmotic contact through a membrane. The glass prism was fitted 
with a pressure transducer for measuring the pressure inside the membrane tube. 

When preparing an experiment, a piece of membrane was fitted in a narrow gap be- 
tween the prism and the cuboid. The prism was then filled to a suitable height with 
a sugar solution and pure water was carefully deposited on top of the sugar solution 
to create a sharp sugar front. Then, the cuboid was filled with water, and the pressure 
transducer was mounted, thereby closing the prism. 

To track the time evolution of the the sugar front inside the prism, we used the re- 
fraction of a laser sheet passing through it. The laser sheet was generated by shining a 
laser beam, generated by a Mcllcs Griot 3.1 mW laser, through a glass rod. When pass- 
ing through the prism, light would deviate depending on the local index of refraction. 
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Figure 4. Experimental results from setup I. Top: Time series of pictures taken in experiment 
5. Time increases from left to right in steps of 30 minutes. See details of the sugar solutions 
used in table [T] Middle, left: Plot of the front position versus time obtained from the images 
above. Bottom, left: Plot of the pressure inside the tube versus time. The (red) dashed line is 
the osmotic pressure of the solution, taken to be the average value of the pressure from t — 2 
h until the end of the experiment. Middle, right: Plots of the sugar front position versus time 
for different sugar concentrations, as indicated in table [T] Bottom, right: Plots of the pressure 
inside the membrane tube for different sugar concentrations. 

The index of refraction varies linearly with sugar concentration and thus by looking at 
the refracted laser sheet projected onto a screen, we were able to reconstruct the con- 
centration profile inside the prism. A camera recorded images of the screen at regular 
intervals to track the moving concentration profile. The procedure used for obtaining the 
concentration profile inside the prism from the images aquired is discussed in appendix 

E2 



3. Experimental results 

3.1. Experimental results, Setup I 

The motion of the sugar front was investigated for solutions of varying sugar concentra- 
tion. An example of a set of data is shown in figure [4] In (A) are the raw images, which 
after processing gives (B) showing the position of the sugar front, function of 

time. The errorbars on Xf are estimated to be ±1 mm, but are too small to be seen. 
Finally, (C) shows the pressure inside the tube as a function of time. At first, a linear 
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1 2 3 4 5 

Mean sugar concentration, c [mM] 1.5±0.3 2.10±0.03 2.4±0.2 4.2±0.7 6.8±0.1 

Osmotic pressure, II [bar] 0.14 ± 0.02 0.15±0.01 0.31±0.03 0.39±0.01 0.68±0.02 

Membrane tube length, L [cm] 28.5 20.8 28.5 28.5 20.6 

Initial front height, I [cm] 4.9 3.7 6.6 6.5 4.8 

Table 1. Data for the experimental runs shown in figure [4] 



motion of the front is observed with a front velocity of ~ 1 cm/h. This is then followed by 
a decrease in the front velocity as the front approaches the end of the tube. The pressure 
is seen to rise rapidly during the first hour before settling to a constant value, indicated 
by a red, dashed line. This constant value is taken to be the osmotic pressure, II, of the 
sugar solution. Looking at (A), one observes that diffusion has the effect of dispersing 
the front slightly as time passes. Below the front, the concentration seems to be uniform 
throughout the cross-section of the tube, and there is no indication of large boundary 
layers forming near the membrane walls. 

Similar experiments with different sugar concentrations were made and a plot of the 
results can be seen on the right in figured] (D) and (E). Qualitatively both the motion of 
the front and of the pressure follow the same pattern as shown on the left. One notices 
that the speed with which the fronts move are related to the mean sugar concentration 
inside the membrane tube, with the high concentration solutions moving faster than 
the low concentration ones. The reason why 2 seems to be moving slower than 1 is 
that experiment 2 was conducted in a slightly shorter membrane tube than 1, thereby 
decreasing the characteristic velocity as we shall see later. 

3.2. Experimental results, Setup II 

Figure \5\ shows the data collected using setup II. At the top, a time series of pictures 
is depicted showing the refracted laser-light projected onto a screen, the time between 
each image being one day. Comparing the upper and lower parts of each picture, one 
generally observes a deflection to the right at the bottom, corresponding to a high sugar 
concentration at the bottom of the prism. In the intermediate region one sees a dip in 
the refracted light, corresponding to a strong concentration gradient. The dip gradually 
flattens while it advances upwards, representing a sugar front which advances while it 
broadens. This process can be seen directly on the left in the middle row, which shows 
the time evolution of the sugar concentration obtained from the images. Starting from 
a steep concentration profile, we see that the front moves forward while it flattens. On 
the left in the bottom row, the time evolution of the concentration gradient is depicted, 
clearly showing a peak which broadens while it moves forward. Finally, on the right, the 
position of the sugar front and the pressure inside the prism is plotted as a function of 
time is shown. The errorbars on Xf are ±1 mm, found as discussed below. 

3.2.1. The effects of diffusion 

To study the effects of diffusion on the dynamics of the sugar front separately, an exper- 
iment was made with setup II, in which the membrane separating the two compartments 
were removed. The experiment was then prepared in the usual way, and the motion of 
the front recorded. The results of this is shown on the right in figure [6] Starting from a 
steep concentration gradient, we observe that the front flattens but otherwise does not 
move much. 

Comparing figures H] and [5] we observe, that while the front moves 2 cm due to osmosis 




Figure 5. Results from setup II. In (A) the raw data images are shown. In (B) the concentration 
profile extracted from (A) is shown. (C) shows the front position extracted from (B) by finding 
the maximum of the concentration gradient, shown in (D). Finally, (E, 1-2) shows the pressure 
inside the prism. The dashed line indicates that the pressure sensor was accidentaly off-line. 
The data shows that the pressure rather quickly (within a few hours) reaches a constant level 
corresponding to the osmotic pressure of the sugar solution (See appendix A 2.2.) 



in 72 hours, it does not move at all in 140 hours due to diffusion. Thus, while diffusion 
has a flattening effect, it plays little role in the forward motion of the front. 

Since the front did not move due to diffusion, the fluctuations in the front position 
seen in figure [6] (D) gives a measure of the uncertainity of a single measurement of the 
front position. Taking the standard deviation of the fluctuations gives an uncertainty of 
±1 mm, shown as errorbars in figure [5] (C). 
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Figure 6. Experimental results from diffusion experiments made with setup II. 



4. Theoretical analysis 

4.1. Front propagation 
Before moving on to a more thorough mathematical analysis we shall, following the 



analysis made by (jEschrich et ali 119721 ) , show that the motion of the sugar front can 
be understood through simple conservation arguments. To that end, let us consider the 
situation in setup I. Let xj denote the position of the sugar-dye front, and let V\ denote 
the volume behind the front and V2 that ahead of it. Taking the tube to be inelastic 
and the fluid inside incompressible, we must have that ^ = — If we let p{x) and 
Po denote the hydrostatic pressure inside and outside the membrane tube and c(x) the 
concentration averaged over the cross-section of the tube at position x, the volume flux 
across a unit area of the membrane is 

J(x) = L p {po - p(x) + RTc(x)) . (4.1) 

Any hydrostatic pressure gradient inside the tube will occur only due to viscous flow, 
but for large tubes and slow flows, this effect is entirely negligible. Thus, since 

L 

J{x')dx' = (4.2) 



we get, that 



p- Po = ^L [ c (x')dx' = RTc (4.3) 
L Jo 
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The rate of change of volume 2 is then 



dt 

Finally, using that 

we get for the front position 



L 



2nrL p I J{x')dx' = -2nrL p RTc(L - x f ) (4.4) 



dV2 2^ X f 

~dt ~ ~dT 



(4.5) 



dxf 2L T) RTc, T . l. T . . , „. 

(L-x f ) = -(L-x f ) (4.6) 



dt r ' 31 t 

which has the solution 



x f (t) =L-(L-/)expl --J (4.7) 

where I = x/(0) is the position of the front at t = 0. This simple result, shows that the 
relaxation time 

(4.8) 



2L p RTc 



for the front propagation depends only on three quantities; the membrane permeability, 
the osmotic pressure of the sugar solution and the ratio of membrane circumference to 
cross-section area. 

The neglect of elastic deformations of the tube is appropriate for the experiments 
made in both setups, as shown in Appendix A. 2, as long as the pressures remain below 
around 1.2 bars. Since we use the large dextran molecules, the osmotic pressures remain 
well below this limit. When sucrose is used, pressures become much larger and elastic 
properties can become important. In Appendix E, we show how the above results would 
change, if clastic properties are taken into account. 

4.2. Derivation of the flow equations 

To formalize the assumptions made above, we will now derive the equations of motion 
for osmotically driven flows, with the geometry of setup I in mind. In appendix [C] we 
shall sec that under certain conditions the equations are also valid in other geometries, 
such as the triangular geometry of setup II. 

The equations of motion for osmotically driven fl ows have be en derived and analyzed 



thoroughly several tim es in the literature (see eg. IWeirl 119811) and h ave been studied 
carefully numerically, ( Thompson fc Holbrook 2003a &T Henton 2002 ). For the sake of 



completeness, we shall include a short derivation of these. 

We consider a tube of length L and radius r, as shown in figure ([7]). The tube has a 
constant cross section of area A — irr 2 and circumference S = 2irr and its walls are made 
of a semipermeable membrane with permeability L p . Inside the tube is a solution of sugar 
and water with concentration c(x). Throughout this paper, we shall study the transient 
dynamics generated by an asymmetrical initial concentration distribution, where the 
sugar is initially localised to one end of the tube with a concentration level Co- The tube 
is surrounded by a water reservoir, modelling the water surrounding the membrane tube 
in setup I. 

We shall assume that L ^> r and that the radial component of the flow velocity 
inside the tube is much smaller than the axial component, as is indeed the case in the 
experiments. With these assumptions, we will model the flow in the spirit of lubrication 
theory and consider only a single average axial velocity component u{x,t). Also, we will 
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assume that the concentration, c, is independent of the radial position, p, an assumption 
that can be verified experimentally in setup II. 

Let us now consider the equation of volume conservation by looking at a small section 
of tube between and Xi. The volume flux into the section due to advection is 



(4.9) 



where the axial flow velocities are taken to be and U{ at an d Xj, respectiv ely. 
The volume flux inwards across the membrane due to osmosis (see eg. Schult d 1980l ) is 



SAxL p {RTc-p). 



(4.10) 



where p is the pressure. For clarity we use the van't Hoff value II = RTc for the osmotic 
pressure, which is only valid for ideal solutions. When comparing with experiments we 
replace this by II = 7c, where 7 must be experimentally determined (see appendix A. 2), 
and thus only use that II is linear in c for small concentrations. Assuming conservation 
of volume, we get that 



A(ui-i -Ui) + SAxL p {RTc - p) = 



(4.11) 



Letting Ax — > and using that the cross-section to perimeter ratio reduces to £, this 
becomes 



r dii 



(4.12) 



For these very slow and slowly varying flows, the time dependence of the Navier-Stokcs 
equation can be neglected and the velocity field is determined by the instantaneous 
pressure gradient though the Poiscuillc or dArcy relation (for a circular tube) 



r dp 
8n dx ' 



(4.13) 



where r\ is the dynamic viscosity of the solution, typically ~ 1.5 x 10~ 3 Pa s in our 
experiments. 

Differentiating equation (|4.12| with respect to x and inserting the result from equation 
(|4.13| we get for the conservation of water that 



RT 



(4.14) 



dc r d 2 u 877 
dx 2L P dx 2 r 2 

The final equation expresses the conservation of sugar advected with veolcity u and 
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Table 2. Values of the parameters M and D in various situations. 



diffusing with molecular diffusivity D 

dc due d 2 c 

The set of equations | 4.14H4.15I are equivalent to those of Thompson and Holbrook 
( Thompson fe Holbrookf 20036 ) except for the fact that we have removed the pressure 



by substitution, and that we do not consider elastic deformations of the tube. 

4.2.1. Non-dimensionalization of the flow equations 

To non-dimcnsionalize equations (|4.14p and (|4.15|) . we introduce the following scaling 

c = cqC, u = uqU, x = LX, t = tor. 

L has been chosen such that the spatial domain is now the unit interval X € [0,1], 
Uo = L/to and cq is the initial concentration level in one end of the tube. Choosing 
further 

r , r 16nL 2 L.„ = D Dr 

t Q = , M= — ^— £ and D = = — , (4.16) 

2L p RTcv r 3 u L 2RTc L 2 L p ' v ; 

and inserting in equations 14. 141 and 14.151 we get the non-dimensional flow equations. 

d 2 U dC 

dx>- MU =dx' (417) 

dC dUC -d 2 C /4<n , 
^ + ^X= D 3X 2 - (418) 
Going back to the original notation 

X — > x, U — » u, C — > c, t — > t 

we finally obtain 

d 2 u dc 

dx~ 2 ~ Mu =dx> ^ 
dc due - d 2 c 

lft + ^- D dx- 2 - ^ 

The parameter M corresponds to the ratio of axial to membrane flow resi stance, which we 

shall re fer to as the Munch number. This is identical to the parameter F in Thompson fe Holbrookl 



(|20036l ). The second parameter D is the ratio of diffusive and advective solute flux. Thus, 
the longer the tube is the less important diffusion becomes and the more important the 
pressure gradient due to viscous effects become. 

Values of the parameters M and D in different situations can be seen in table [2] The 
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Quantity 

Radius [/im] 

Length [mm] 

Flow velocity [mh _1 ] 

Elastic Modulus [MPa] 

Permeability [1CP 11 ms~ Pa~ 

Sucrose concentration [M] 



Magnitude Reference. 
4.5 (Fava bean), 4 (Winter squash), 6-25 f, *, | 

0.09 (Fava bean), 0.1-3 f, II 

0.5-1, 0.2-2 f, || 

17, 5.6-7.4 (Ash) •, ( 

5,1.1 (Zitella translucence) •, ) 

0.3-0.9 * 



Table 3. Characteristic p roper ties of phloem siev e el ements. References: flK noblau ch fc van Bel 
' 19981) . * iTaiz fc Zeigerl (l2002h . || iNobell l| 19991 ). • iThompson fc Holbrookl (|2003at ). { iNiklas 



1992), ) Esch rich et al\ (|1972 



typical magnitude of the parameters M and D in plants are found from the values also 
given in table [2j 

r = 10 fixa, i] = 1.5 x 1(T 3 Pas, u = 2 mh" 1 , L p = 2 x 10 _11 m(Pas)" 1 . 

We observe, that M and D are small in both experiments, and that for short distance 
transport in plants this is also the case. However, over length-scales comparable to a 
branch (L = 1 m) or a small tree (L = 10 m) M is large, so in this case the pressure 
gradient is not negligible. 

When deriving the equations for osmotically driven flows, we have assumed that the 
concentration inside the tube was a function of x and t only. However, the real con- 
centration inside the tube will also depend on the radial position p in the form of a 
conce ntration boun dary layer near the membrane, in the literature called an unstirred 
layer (|Pedlevlll983h Close to the membrane, the concentration c m is lowered compared 



to the bulk value, q,, because sugar is advected away from the membrane by the influx 
of water. This, in turn, results in a lower influx of water, ultimately causing the axial 
flow inside the tube to be slower than expected. In our experiments we see no signs of 
such boundary layers and apparently their width and the effects on the bulk flow are 
very small. 



5. Solution of the flow-equations 

The equations governing the time evolution of a sharp sugar front has been known for 
the closed and semi-closed tube since the work of Eschrich et al. . However, their solutions 
yield only the position of the front, and not the concentration profiles in front of and 
behind the concentration front. Analytic solutions giving the time evolution of the entire 
concentrati on profile has been fo und for the c losed tube for M = D = by G. J. Weir 
(|Weirlll98ll ) and by H. L. Frisch (|Frischlll976T ) for the the semi-closed tube for M = 0, 
D ^ 0. In both cases the authors have started from piecewise constant Heaviside-like 
initial concentration profile. To extend this work, we shall present a method for solving 
the equations of motions analytically using Riemann's method of characteristics. For an 
arbitrary initial condition, this method will generally yield an implicit solution. Only 
in special cases will it yield a closed formula for the solution. The method works for 
M = D = and for closed and semi-closed tube geometries. 

For arbitrary values of M and D, we cannot solve the equations of motion analytically 
and thus have to use nu merical methods. This topic has been the focus of much w ork both 
in the steady-state case ( Thompson fc Holbrookf2003a ) and in the transient case ( Hentonl 
20021) . However, no formulation capable of handling all different boundary conditions has 



so far been presented. Therefore, we show that using Green's functions, the equations 
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of motion can be transformed into a single integro-differential equation, which can be 
solved using standard numerical methods. 



5.1. Results for small Munch number 
In the limit M = D — the equations become 

d 2 u dc 
dx 2 dx ' 



dc due 



0. 



dt dx 

By integrating equation equation 15. II with respect to x, we get that 

du , . 

Tx- C + F{t) - 

If we choose u(0) — u(l) = 0, F(t) becomes 



F(t) 



c dx = —c(t), 



Using 15.31 in equation 15.21 gives 



d /du u (® u 
dx \ dt \ dx 



(5.1) 
(5.2) 

(5.3) 
(5.4) 
(5.5) 



Integrating with respect to x and using the boundary conditions on u, this becomes 



du du 
dt dx 



-cu. 



(5.6) 



Equation [STB] is a damped Burgers equation (see eg. lGurbatov et aZ.lll991r ). which can be 
solved using Riemann's method of characteristics. The characteristic equations are 



d.c 
dt 



-cu 



(5.7) 
(5.8) 



Equation 15.71 has the solution 

u = uo(£) exp(— ct), (5-9) 
where the paramctrization £(x, t) of the initial velocity has to be found from 



x = £ + -u (0 (1 - exp(-ct)) 
c 



(5.10) 



where ^ = x at t = 0. 



5.1.1. Solution for piecewise constant initial concentration 

To be able compare our method to the results obtained by (|Weirill98lh . we will use a 
Heaviside step function as initial condition on c 



c(x,t = 0) = aH{\ - x) 



d for < x < A, 
for A < x < 1, 



(5.11) 
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Equation 15.31 now enables us to find the initial condition on the velocity 

u(x,t = 0) = (c(a/,0) -c) dx' = / (c(x', 0) - Ac 7 ) dx' (5.12) 
Jo Jo 



(ci — c)x for < x < A, 
c(l — x) for A < x < 1, 



(5.13) 



From eguation l5.131 we have that 



«o({) = | ( ?-f °" e - A ' (5.14) 

W \c(l-0 for A<C<1. 

Then, solving for £(x,t) in equation [STTU] gives 

C( ' ^l^r*^ ^ ae / a (5 - 15) 

where the intervals I± and li are defined by 

h = [0,1- (1 - A)exp(-ct)], (5.16) 

J 2 = [l-(l-A)exp(-ci),l]. (5.17) 
Finally, u(x,t) is calculated from equation 15.91 

*(l-A)(l-e*p(- Et ))' (518) 

c(l — x), for x £ ^2 



which is equivalent to the result obtained by (|Weirlll98lh . The solution is plotted in figure 



[SI top. We can now calculate the instantaneous sugar front position Xf and velocity Uf 
using the right boundary of I\ from equation 15. 161 

x/(i) = l-(l-A)exp(-cf), (5.19) 

Uf (t) = — j- = c(l - A)exp(-ci). (5.20) 
Similarly, c(x, <) is given by 

1 — (1 — A) cxp(— ct) 
Going back to dimensional variables, equations (|5 . 19[) and (|5.20[) become 

X f (t) =L- (Z-Z)exp (-- ) and (5.22) 



u f (t) = -exp(--\, (5.23) 



r \ r 



where L is the length of the membrane tube, I is the initial front position and the decay- 
time r is accordance with the simple argument leading to equation 
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f 

Figure 8. Top: Plot of the analytical solution for a piecewise constant initial concentration. 
A = 0.1, ci — 1 and c = 0.1. Bottom: Plot of the analytical solution for a piecewise linear initial 
concentration. Ai = 0.05, A2 = 0.15, a = 1 and c = 0.1. Time increases from black to gray in 
steps of one unit of time. 



5.1.2. Solution for piecewise linear initial concentration 

As initial condition on c, we will use the piecewise linear concentration profile given 

by 

!ci for < x < Ai, 

c/xSr for Al - x - A2 ' ( 5 - 24 ) 
for A 2 < x < 1, 

Using [5731 yields the initial velocity 

!{ci — c)x for < x < Ai, 

A lX 2 + B x x + d for Ai<a;<A 2 , (5.25) 
c(l - x) for X 2 < x < 1, 

where c = ci Xl + X ' 2 , and the constants are given by 

(5.26) 



Finding from equation 15. 101 now gives, that 



l+i(l-A)(l-cx P (-et) 



for a; € Ii, 



^{x, t) = I A 2 H + B 2 H + C 2 for x G h, (5-27) 



x — l+cxp( — ct) 
cxp( — ci) 



for X € I3, 
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where 

A B 

A 2 = -^(l - exp(-ct)), B 2 = 1 + — (1 - exp(-ci)), 
c c 



C 2 = ^(l-exp(-a)), 
(5.28) 



and 



fa 



-B 2 + ^B 2 -4A 2 (C 2 -x) 
2Ao 



where the plus solution has been chosen to ensure, that £ — > x as < — ► 0. Finally, 



/a 



0, Ai + — (a - c)(l - exp(-ct)) 
c 

Ai 



Ai + — (cj — c)(l — exp(— ct)), 1 — (A 2 — 1) cxp(— ct) 
c 

h= [1 - (Aa - 1) exp (-ct) , 1] . 
Plugging into 15.91 gives u(x,t) as 



(cj — c) exp( — c)a: 
l+i(l-A)(l-cxp(-ct) 



for a; G 7xi 

*) = (Ai^f + SiCl + Ci) exp(-a) for a? e la, 

t c(l — a;) for a; e ^3, 

as shown in figure [5] along with c found from equation [57 



(5.29) 

(5.30) 

(5.31) 
(5.32) 

(5.33) 



5.2. Results for large Munch number 

In the limit of large M > 1 we cannot neglect the pressure gradient along the channel 
and this term dominates the advective term in (|4.14[) , i.e. the second derivative in u. 
Thus 

dc 
Ox 



-Mu 



(5.34) 



dc dcu 
dt dx 

jiving the nonlinear diffusion equation 

dc ^ d dc 
dt dx dx 



D 



d 2 c 
dx 2 



D 



d 2 c 
dx 2 



(5.35) 



(5.36) 



If we neglect molecular diffusion, which is true as long as Mc ^> D « 10 , the resulting 
universal nonlinear diffusion equation can be written 



^ C — M ® 
dt dx 



dc 
dx 



(5.37) 



which belongs to a class of equations which have been studied e.g. in the context of 
intense thermal w aves by Zeldovich et al. and flow through porous media by Barenblatt 
( Barenblatt 19961 ) in the 50'ies. The Munch number M can be removed by rescaling the 
time according to T = Mt, so in this limit we get very slow motion with a time scale 
growing linearly with M. The equation (|5.37|) admits scaling solutions of the form 



c(x,t) = (Mt) a $(0 with f = x(Mtf 



(5.38) 



18 



K. H. Jensen, E. Rio, R. Hansen, C. Clanet and T. Bohr 



as long as a + 2/3 + 1 = 0. The total amount of sugar is, however, conserved. In our 
rescaled units 

c(x) dx = X (5.39) 



where, as before, A is the fraction of the tube initially containing the sugar. We can only 
hope to find a scaling solution in the intermediate time-regime, where the precise initial 
condition has been forgotten, but the far end (x = 1) is not yet felt. Thus we can replace 
the integral (|5.39[) with 



cdx = X (5.40) 



o 



which implies that a = (3 = — 1/3 and 



c(x,t) = (Mt)- 1/3 <i>(£) with £ = — ^ (5.41) 

Inserting this form into (|5.37[) . we obtain the differential equation for <i> 

1 d 2 <f> 2 1 d(£$) 



2 d£ 2 3 d£ 

which can be integrated once to 



(5.42) 



d$ 1 

$ h -£$ = const (5.43) 

d£ 3 

Due to the boundary condition dc/dx = in the origin, the constant has to vanish and 
we find the solution 

m = \(b 2 -e) (5.44) 

which is valid only for £ less that the constant b. For £ > 6, $ is identically 0. The fact 
that the solution - in contrast to the linear diffusion equation - has compact support, is 
an in teresting characteristic of a large class of nonlinear diffusion equations (jBarenblatt 



19961) . The value of b is determined by the conservation integral (|5.40[) giving / °° $ d£ 
1, and thus 6= (9A) 1 / 3 . 

The final solution thus has the form 

c(M) = /^(M0) 2 -* 2 ) for x < x f (t) = (9AMt)V3 
1 for x > Xf(t) 

which shows that the sugar front moves as Xf(t) ~ t 1 ^ 3 and the concentration at the 
origin decays as c(0, t) ~ i -1 / 3 . To check the validity of this solution, also when the initial 
condition has support in a finite region near the origin, we plot (cL)~ 2 ^ 3 (Mt) 1 ^ 3 c(x, t) 
against £ = x^cLMt)" 1 ^ 3 in figure [H C. 

5.3. Numerical methods for nonzero M and D 
For nonzero values of M and D, the equations of motion, 



and 

dt ' dx dx 



dc . dcu = D Pc (g 4?) 
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M=100 




Figure 9. (A) Numerical simulation of equation (|5.37p compared with (Bj the scaling solution 
(|5.45j) and (C) (|5.44|) . The intial condition has the form c(x, 0) = 1 — — ( *-xj where A = 0.1 

and e = 2 x 1CP 2 and the curves are equidistant in time. When A controlling the size of the 
region of nonzero initial sugar concentration becomes larger, a more accurate scaling solution is 
found by letting t — > t + to and treating to as an unknown parameter. In (C), we have omitted 
the first curve (the initial condition). 



cannot be solved analytically. However, they can be written as a single integro-differential 
equation, which is straightforward to solve on a computer. If we choose a set of linear 
boundary conditions, -B^M = asj, for equation 15. 46i the solution can be written as 

G(x,0^dt + u 2 . (5.48) 



o 



Here, G(x, £) is the Green's function for the differential operator — M with boundary 
conditions B x [u] = and u 2 fulfills the homogeneous version of 15.461 with B x [u] = a; . 
Plugging this into equation 15.471 yields 



For the closed tube, ie. for the boundary conditions u(0, t) = u(l, t) = 0, G(x, £) is given 
by 

l-ffsff smh(a(l-0) for g > x, 

and ui = 0. To increase numerical accuracy, it is convenient to transform equation 15.491 
by defining 

g = c-a (5.51) 
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Figure 10. Results of numerical simulation of equation 15.491 using the boundary conditions 
u(0,t) — u(l,i) — for different values of M. D is kept constant at 10~ J . The initial condition 
was c(x, 0) = 1 r where A = 0.2 and e = 2 x 10 2 

l + cxp(^ — J 



and choosing /(0) 
we get that 



/(l) = such that f(x) = f (c — c) <i£. Inserting in equation 15.491 



df 



l dK(x 2 e L 



df 

dx 



(5.52) 



where 



dK{x,t) 



~a binh ^} x ^ sinhaf for £ < x, 

sinn a ^ ■» j 

-«^sinh(a(l-£)) for f > x. 



(5.53) 



To solve eauation l5.52l we used Matlab's built-in time solver ode23t which is based on 
an explicit Rungc-Kutta formula along with standard second order schemes for the first 
and second order derivatives. For the spatial integration, the trapezoidal rule was used 
(jPressI 120011 ) . The numerical code can be found in appendix [Dl Results of a numerical 
simulation for different values of M is shown in figure 15.31 
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Figure 11. (A): Experimental (black dots) and fits to equation l5,22l for the relative front 
position vs. time. (B): Semi-logarithmic version of (A). 
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Figure 12. Our e xperimentally obtain ed values of to plotted together with the results found 
by Eschrich et al. (|Eschrich et fflZ.lll972f ). Data points marked with an a represents results from 
closed tube experiments and points marked with a b represents results from semi-closed exper- 
iments take from the original paper, figures 8 and 9. 



6. Comparison between theory and experiment 

In section [3J we have presented experiments demonstrating the movement of a sugar 
solution inside a membrane tube surrounded by a reservoir of water. We now wish to 
consider whether the theory is in agreement with the experimental results. 

6.1. Setup I 

The plot in figure [Til shows the relative front position, *f , plotted against time for 
five different experiments conducted with setup I. The numbers 1 — 5 indicates the sugar 
concentrations used, cf. table [T] One clearly sees, that the relative front position ap- 
proaches zero faster for high concentrations than for low. Typical values of M and D 




are M ~ 10 -8 and D ~ 10~ 5 , so it is reasonable to assume that we are in the domain 
where the analytical solution for M = D = is valid. To test the result from equa- 
tion 15. f 91 against the experimental data, the plot in figure [TT] shows the logarithm of 
the relative front position plotted against time. For long stretches of time the curves are 
seen to approximately follow straight lines in good, qualitative agreement with theory. 
The red dashed lines are fits to equation (|5.19|) . and we interpret the slopes as — the 
different values plotted in figure [T2l against the theoretical values. The theoretically and 
experimentally obtained values of to are in good quantitative agreement, within 10-30%. 
Generally, theory predicts somewhat smaller values of to than observed, implying that 
the observed motion of the sugar front is a little slower than expected from the pressure- 
flow hypothesis. Nevertheless, as can be seen in figure [T2l these results are a considerable 
improvement to the previous results obtained by Eschrich et al. as we find much better 
agreement between experiment and theory. 



6.2. Setup II 

The plot in figure [13] shows the relative front position, L *f , plotted against time for 
the experiment conducted with setup I. On the semi-logarithmic plot, the curves are 
seen to follow straight lines in good, qualitative agreement with the simple theory for 
M = D = 0. As can be seen in figure[l2l we also found very good quantitative agreement 
between the experiment and theory for setup II. 

To test how well the motion of the sugar front observed in the experiments with setup 
II was reproduced by our model, we solved the equations of motion numerically starting 
with the initial conditions from figure El For M = D = 0, The results are shown as (red) 
curves in figure[14](B). While the front positions are reproduced relatively well, the shape 
of the front is not, so diffusion must play a role. This can be seen in figure [Til (C) which 
shows the result of a simulation with M = 10~ 9 , D = 6.9 x 10 _1 m 2 s _1 . Clearly, the 
model which includes diffusion reproduces the experimental data significantly better. 

To study the shape of the front in greater detail, consider the plots on the right in 
figure [T5J Here the gradient of the concentration curves on the left in figure [TJ] arc shown. 
In (A) we clearly see a peak moving from left to right while it gradually broadens and 
flattens. In (B) we also see the peak advancing, but the flattening and broadening is 
much less pronounced. In (C) we see that the model which includes diffusion reproduces 
the gradual broadening and flattening of the front very well. 
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Figure 14. Results from setup II showing the experimental data (top row) and the numerical 
model for M = D = (middle row) and for M = 1CT 9 , D = 6.9 x Kn u m 2 s _1 (bottom row). 

7. Summary 

In this paper we have studied osmotically driven, transient pipe flows. The flows are 
generated by concentration differences of sugars in closed tubes, fully or partly enclosed 
by semi-permeable membranes surrounded by pure water. The flows are initiated by a 
large concentration in one end of the tube and we study the approach to equilibrium, 
where the sugar is distributed evenly within the tube. Experimentally, we have used two 
configurations: the first is an updated version of the setup of Eschrich et al. where the 
flow takes place in a dialysis-tube and the sugar is followed by introducing a dye. The 
advantage is the relatively rapid motion, due to the large surface area. The disadvantage 
is that the sugar concentration cannot be inferred accurately by this method and for 
this reason we have introduced our second setup, where the sugar concentration can be 
followed directly by refraction measurements. On the theoretical side, we first re-dcrive 
the governing flow equations and introduce the dimcnsionless Munch number M. We 
then show that analytical solutions can be obtained in the two important limits of very 
large and very small M. In the general case we show how numerical methods based on 
Green's functions are very effective. Finally, we compare theory and experiment with 
very good agreement. In particular the results or the velocity of the front (as proposed 
by Eschrich et al.) can be verified rather accurately. 

It is a pleasure to thank Francois Charru, Marie-Alice Goudeau-Boudeville, Herve 
Cochard, Pierre Cruiziat, Alexander Schulz, N. Michelle Hollbrook and Vakhtang Putkaradzc 
for many useful discussions. Much appreciated technical assistance was provided by Erik 
Hansen. 



Appendix A. Setup and Methods 

A.l. Materials 

A. 1.1. Chemicals 

The sugar used was a dextran (Sigma- Aldrigde, type D4624) with an average molecular 
weight of 17.5 kDa. The dye used was a red fruit dye (Flachsmann Scandinavia, R0d 
Frugtfarve, type 123000) consisting of an aquous mixture of the foo d additives E-124 and 
E-13 1 with molecular weights of 539 Da and 1159 Da respectively (jPubChem-Database 
Even though the molecular weights arc below the MWCO of the membrane, the 



red dye were not observed to leak through the membrane. This however, was observed 
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Figure 15. Left: van't Hoff relation for 175kDa Dextran. Right: Elastic properties of the 

membrane tube 



when using another type of dye, Methylene blue, which has a mulecular weight of 320 
Da. 

A. 1.2. Membrane 

The membrane used in setup I was a semipermeable dialysis membrane tube (Spec- 
tra/Por Biotech cellulose ester dialysis membrane) with a radius of 5 mm, thickness 60 
/xm and a MWCO of 3.5 kDa. The membrane used in setup II had identical specifications 
except that it had a radius of 3 cm allowing it to cover the interface between the prism 
and the water reservoir after beeing cut in half. 

A. 2. Elastic and osmotic properties of the materials used 

A. 2.1. Elastic properties of the membrane tube 

Figure [12] (right) shows the relation between internal pressure and radius for the mem- 
brane tube. For pressures less than 1.2 bar, a linear relation betwe en the relat ive radial 
increase and internal pressure was found. Linear elasticity theory ( Love! 1944 ) predicts 
that for a thin-walled cylindrical tube 

r = r„ + f, (Al) 

where ro is the equilibrium radius and d the thickness of the tube, p is pressure and E 
is Young's modulus of the membrane tube. From this, E, was found to be 

E = 0.66 ± 0.01 GPa (A 2) 

A. 2. 2. Osmotic strength of Dextran 

Figure [15] (left) shows the relation between dextran concentration and osmotic pressure 
found from the experiments shown in figured] A linear fit gives 

n = (0.1 ±0.01 bar mM _1 )c (A3) 

where LI has un its of bar, and c is measured in mM. This is in good agreement with 
values given by ( Jonsson 19861 ) 



Appendix B. Tracking the sugar front position 

B.l. Setup I 

After running an experiment with setup I, the raw data we had acquired consisted of a 
series of pictures as shown in figure |U 
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Figure 16. Data processesing in setup I. 

(A) : Raw RGB data image showing the sugar solution mixed with red dye inside the membrane 
tube. 

(B) : The image Irom (A) has been cropped and filtered. 

(C) : The solid black curve is the intensity of color in (B) taken along the white vertical line. 
The red curve is the absolute value of the gradient of the intensity. The front position is taken 
to be where this curve has its maximum value (at approximately 200 pixels). The black dot in 
(B) is the position of the front found in this manner. The peak in the intensity gradient near 
920 pixels is due to the membrane fitting, and was ignored. 

To track the position of the sugar-dye front, the images was treated as shown in figure 
[T6l First, the image was imported into Matlab as (A), and then cropped to shown only 
the membrane tube, (B). Simultaneously, it was filtered to give the highest contrast for 
obtaining a well-defined front position. To find the front position, a vertical line running 
along the center of the membrane tube was picked out, shown as a white line in (B). 
Along this line, the color intensity was found, shown as the black curve in (C). Finally, 
the gradient of the color intensity was found - shown in (C) as the red curve - and the 
front position was defined to be the the position of the maximum in intensity gradient. 

To justify the use of red dye as the tracking medium, we took closeup images of the 
sugar front as shown in figure [T71 It is clearly seen, that the dye moves with the point at 
which the concentration gradient is largest. Thus, we conclude that the dye travels along 
with the sugar. 

B.2. Setup II 

A camera recorded images of the screen at regular intervals, as shown in figure El The 
deflection at the bottom of the image corresponds to a high sugar concentration inside 
the lower part of the prism, and the vertical deflection is due to a strong gradient in 
index of refraction near the sugar front. 

As the beam passes through the prism, it gets deflected due to variations in index of 
refraction of the fluid inside the prism relative to the surrounding air. To determine the 
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Figure 18. Data processing in Setup II. (A): Raw data image. (B): Filtered data image. (C): 
Index of refraction inside the prism giving the red curve in (B) cf. equation lB5l and IbTuI 



horizontal deflection of a light beam passing through the prism, we consider the situation 
sketched in figure [TU The deflection Ai along the y axis is given by 

Ai=Gtan6i (Bl) 

where G is the orthogonal distance from the prism to the screen, and 0i is the deflection 
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Figure 19. Deflection of the laser beam as i passes through the prism. 



angle as shown in the figure. To find 0i we notice that 

n -6i = 



(B2) 



where ^ is t he prism angle. The deflection angle, ri, is given by the Sncll-Descartes law 
(|Hechtll2002j ) 

n sin — sin r\ , (B3) 

where n is the index of refraction of the fluid inside the prism. In the experiments, 0i is 
typically small so 



nsin^/ = sinri = sin(^E' + Oi) ~ sin^ + Q\ cos 'I', 
so the horizontal deflection 

Ai ~ G6i =G(n- l)tan*. 



(B4) 
(B5) 



To determine the vertical deflection consider the situation sketched in figure [T§I (B). 
According to Format's principle, light travels along the path that can be traversed in the 
least possible time. Consequently, as the light passes through the prism, it travels in a 
circular arc with a radius of curvature, 1Z, given by 



K n ' 



(B6) 



where N is the normal to the trajectory of the light beam lLandau fc Lifshita (|1984l ). The 
vertical deflection is given by 

A 2 = Gtane 2 (B 7) 
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where the angle Q2 is given by 

sin 2 = n sin r 2 (B 8) 

Since r 2 is small, it is clear from the figure that sin r 2 ~ ^ , so from equation (|B 6p 

•9 9 dn 

SUlr ^TZ = nd^- (B9) 
In the limit where r 2 and 2 are small we get for the vertical deflection 

A 2 =G.g^, (BIO) 

where we have assumed that Vn is constant as the beam traverses the prism. Having 
obtained Ai and A 2 we now have to deduce n(x) from the projected image. We will do 
this by assuming that n(x) has a generic, sigmoid shape 

n(x) = no + ni(l-— x _ t ) (Bll) 

where the constants no and n\ controls the magnitude of n and / and e controls the 
position and steepness of the front. A plot of this function can be seen in figure (j!8l C) 
The procedure for obtaining n(x) was as follows. First, the raw data image was loaded 
into Matlab as shown in figure [T8l (A). Then, the image was filtered to show only 
regions of high light intensity, shown as the black dots in (B). Then, a guess of the the 
form IB 1 T1 was made, and the deflections Ai and A 2 was calculated from equations IB 51 
and IB 101 shown as the red curve in (B). Finally, an optimization of the parameters was 
made using Matlab's fminsearch engine, thereby giving the n{x) of the form IB 111 best 
able to reproduce the image seen in (A). Generally the assumption that n was of the 
form IB lH gave very good fits, as can be seen in figure IT51 (B). 



Appendix C. Generalization of the equations of motion to 
non-cylindrical geometries 

When deriving equations 14. 19114.201 we have assumed that our system consisted of a 
cylindrical tube with semipermeable walls. The assumption of a cylindrical tube, however, 
have only been used in equation (|4.12[) where we assumed that the cross-section area to 
perimeter ratio was 

H (C1) 

and in equation (|4.13|) where we assumed the axial resistance in Stokes flow to be inversely 
proportional to the cross-section area 

2tt _ 8 
A ~ r 

These two factors are of purely geometrical nature and appear in M, D and <o a s 

Moc---! = ^, (C3) 

D^ T - (C4) 

and 

i cx^. (C5) 



(C2) 
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Therefore, as long as the assumption of a ID flow velocity and concentration holds 
inside the tube the equations of motion can be extended to include other geometries, e.g. 
triangular tubes as used in setup II, by replacing the geometric factors in M, D and to 
as discussed above. 

Finding the cross-section area to perimeter ratio is trivial, and the expression for the 
axial resistance in Stokes flow generally has the form 

lAdp fnf(\ 
u= (C6) 

77 a ox 

where ^ is a purely geometric f actor, which for a cylind rical tube is . For various pipe 
cross-sections Mortcnsen et al. (jMortensen et aLll2005l) has found a as a function of the 
dimcnsionless compactness 

02 

C=-. (C7) 
C.l. Setup II 

To extend the equations of motion to setup II, consider the following. For a isosceles 
right triangle with two sides of length s and one of length \^2s, we get that 



2(2S+ /^12 + 8V2- 



C* = ^ ^ — ^ = 12 + 8V2 (C8) 

s z 

Mortcnsen et. al. showed that for pipes with triangular cross-sections 

« = §C+^ (09) 



so m our case 



300 , 200V2 + 40V3^ 38 _ 36 _ 



Also, 



17 17 17 



.4 



s 4 + 2V2" 

Plugging into the expressions for Al, D and to we get 



(Cll) 



= 38.36(8 + 4^L% = ^V^L, (c u) 

D n = (2 + V2)Ds Ds 

(A + 2V2)2RTc L 2 L p 2RTc L 2 L p 

t u_ (2 + V2)s _ s 

(A + 2V2) L p RTc Q 2L p RTc ' v ; 

The extra factors of 2 + y/2 comes from the fact that the membrane only covers one wall 
of length a, thereby scaling L p down to 2 ^^ - 



and 



Appendix D. Numerical code 

1111 25 Feb 2008 V/XL 

1111 Numerical code for solving the equations of motion '/,'/,'/,'/, 
1111 for osmotically driven flows in a closed tube VLVL 
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function numsolver 

'L'L'L'L Choose M and D 'L'L'L'L 
D = le-8; 
M = le-8; 
a=sqrt(M) ; 

VLVL Initialize x and t VLVL 

Nx = 500; 7. 

x = linspace(0, 1 ,Nx) ; 

Nt = 100; 

tspan = linspace(0,40,Nt ) ; 

'L'L'L'L Initial condition TLTL 
epsilon=2e-2 ; 
x0=0 . 2 ; 

c0=l-l./(l+exp(-(x-x0)/epsilon )) ; 
cbar=trapz (x,c0) ; 
for i=2:Nx 

f0(i) = trapz(x(l:i) ,c0(l:i) .cbar) ; 

end 

'L'L'L'L Run solver 'L'L'L'L 

[t,f] = ode23t (@f deriv, tspan, fO, [] ,D, a, cbar ,x,Nx) ; 

'L'L'L'L Plot results TLTL 
TL Calculate c °L'L 
for i=l:Nt 

c (i , : )=gradient (f (i , : ) ,x)+cbar ; 

end 

7.7. Plot c 7.7. 
figure (1) 
surf (c) 

TI'I'L Function fderiv for use in ode23t T/'/°L 
function df = fderiv(t,f ,D,a,cbar,x,Nx) ; 
P = zeros(Nx, 1) ; 
f(l) = 0; 
f(Nx) = 0; 
temp = f ' ; 
f = temp; 
7.7. Calculate P 7.7. 
for i=2:Nx-2 

P(i) = f (i) . (a/sinh(a)) . (trapz (x(l : i) , . . . 

sinh(a. (l.x(i))) . .f (l:i) 

sinh(a. . x(l : i) ) )+trapz (x(i+l : end) , . . . 

sinh(a. . (1 . x(i+l : end) ) ) . .f (i+l:end) 

sinh(a.x(i)))) ; 
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end 

°/o°/ Impose boundary conditions on P %7o 
P(l) = 0; 
P(Nx) = 0; 

P(Nx-l) = P(Nx-2)/2; 
P=P> ; 

laplacef = 4.del2(f,x); 
laplacef (1) = 0; 
laplacef (Nx) = 0; 
7,°/, Update f 11 

df =D*laplacef -P* . (cbar+gradient (f , x) ) ; 
°/ °/o Impose boundary conditions on f '/,'/, 
df(l) = 0; 
df(Nx) = 0; 



Appendix E. Front propagation in an elastic tube 

In this section we shall investigate the effect of expansion of the tube. First we simply 
assume that water can enter the tube and make it expand without an important increase 
of pressure. The total volume of the tube is V = irr 2 L and the volume of the part with 
sugar of concentration co is V — irr 2 x. In the simplest case, we assume that the flow 
inward only up in he sugar interval [0, a;]. The rest of the tube just expands without any 
restoring force (E = 0) 

dV „ T dr 2 dx „ dr 

= 2?rLr — = irr 2 h 2nrx — (El) 

dt dt dt dt y ' 

Thus 

(L - x)r 2 = (L - x )r 2 (E2) 

The inflow is 



I = = 2nr [ L w RTc{x')dx' = 2irrL w RTcx 

dt Jg 

But c changes so that cxr 2 = c^x^r^ and thus 



(E3) 



dV dr 2TTL w RTc x rl 

-dT = 2 * Lr dt = r (E4) 

so 

dr L w RTcoXorQ 
dt, ~ Dr 2 

Note that cq os not the average over the whole tube. This would be c"o = We have 

earlier introduced the time scale 



(E5) 



(E6) 



Using instead 

letting t — sf and r = yro we get 



2L w RTc 



(E7) 



2L w RTcq 
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with solution 

/3s \ 1/3 

y(s) = (- + l) (E9) 

Now x can be found from (IE 2[) : 

(L - x)y 2 = (L — x ) (E 10) 

so 

3s \~ 2/3 f J ^ ~ 2/3 



x = L-{L-x )[— + l\ = L-(L-x ){— + l) (Ell) 



Now the front velocity is 



v{t ) = x>(t) = ^—^^ + i y t/3 (E12) 
and initially it takes the value 

Vo=x i(Q) = L^± (E13) 

T 

We now assume that the influx of water creates an increased pressure, and that, conse- 
quently water can flow out of the sugarless region. We assume that the elastic properties 
are governed by 

p(r) = d Q E r —p± (E 14) 

where do is the initial width of the membrane, E is Young's modulus and we have assumed 
incompressibility (i.e. Poisson's ratio v = 1/2) and, so dr = d^r^. Now the equation for 
the current is 

dV f L 

I = — = 2nrL w / (RTc(x') — p(x')) dx' = 2nrL w (RTcx - pL) (E15) 
at J 

where we have assumed (as in our experiments) that the pressure is constant along the 
tube. Thus 

dV ( L(r-r )\ 

— = 2irrL w iRTcx- doE v °M (E 16) 

Again, cxr 2 = coXoTq so 

— = 2nLr- = 2nrL w d E—^- j (E 17) 

so 

dr _ L w RTc r 2 _ L w d E(r - r ) 

dt r 2 r 2 

Using again r = roy and t = fs we get the dimensionless equation 

dy = 1 f l-B{y-l) 
ds 2\ y 2 

where 

B= 2fdoEL w = do^_ 
To ro RTco 

In this case we see that the radius saturates for long times to the value: 

1 -)- B 

rf = r y f = ro— — (E21) 



(E18) 
(E19) 
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and B can be expressed as 

B= — - — (E22) 
Vf - 1 

This equation can be solved for s{y) by inverting: 

ds v 2 ,„ 

— = 2 E 23 

dy l + B-By V ' 

which can be integrated (with the initial condition s(l) = 0) to 

<V) = ^3 (-( B (y - !)( 2 + B ( 3 + V)) + 2 ( X + B f lo g(l +B- By)) (E 24) 

valid as long as y < yf. For small B this is is approximately: 

s(y) = l/3(-l + y 3 ) + 1/12(1 - Ay 3 + 3y 4 )B + o(B 2 ) (E 25) 

and thus agrees with the result of last section, i. e. the result for E — > 0. For larger B 
the part of the integral close to 

To find the front velocity Vf = x we look at the volume V2 = Trr 2 (L — x) above the 
sugar front. By definition, the sugar concentration is zero there, so the water flux through 
the membrane is only due to the pressure difference: 

^ = 2nr(L - x)— - nr 2 — = -p(L - x)L w 2nr (E 26) 

dt x ' dt dt K ' v ' 



so 



or 



dx _ dlog(L - x) _ 2 fdr , ^^(r-roT. (E2?) 



L — x dt dt r \dt 

d\og(L~x) 2(dy _ 2t Jy-l) 

3 = — [ ~r + Tr o L w 2d E — 

ds y \ds y^ 

V \ds y 



-A^ + B^) (E28) 



from which we can find x(t) from the solution for r(t), although we cannot write it 
explicitly since the explicit form of r(t) is not known (only the inverse (|E 24|) V For B = 
we again recover the fiappy limit (E=0), where L — x ~ y 2 - 
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Using Stokes flow between eccentric counter-rotating cylinders as a prototype for bounded 
nearly parallel lubrication flow, we investigate the effect of a slender recirculation region 
within the flow field on cross-stream heat or mass transport in the important limit of high 
Peclet number Pe where the 'enhancement' over pure conduction heat transfer without 
recirculation is most pronounced. The steady enhancement is estimated with a matched 
asymptotic expansion to resolve the diffusive boundary layers at the separatrices which 
bound the recirculation region. The enhancement over pure conduction is shown to vary 
as e 1 / 2 at infinite Pe, where e 1 / 2 is the characteristic width of the recirculation region. 

I Icy 

The enhancement decays from this asymptote as Pe ' . 



1. Introduction 

The use of integral equations to solve 'exterior' problems in linear acoustics, i.e. to solve 
the Hclmholtz equation (V 2 -I- k 2 )<fi — outside a surface S given that </> satisfies certain 
boundary conditions on S, is very common. A good description is provided by ?. fntegral 
equations have also been used to solve the two-dimensional Helmholtz equation that 
arises in water-wave problems where there is a constant depth variation. The problem of 
wave oscillations in arbitrarily shaped harbours using such techniques has been examined 
(see for example ??). 

In a recent paper ? have shown how radiation and scattering problems for vertical 
circular cylinders placed on the centreline of a channel of finite water depth can be 
solved efficiently using the multipole method devised originally by ?. This method was 
also used by ? to prove the existence of trapped modes in the vicinity of such a cylinder 
at a discrete wavenumber k < n/2d where 2d is the channel width. 

Many water- wave/body interaction problems in which the body is a vertical cylinder 
with constant cross-section can be simplified by factoring out the depth dependence. 
Thus if the boundary conditions are homogeneous we can write the velocity potential 
4>(x,y,z,t) = Re{(j)(x, y) coshfc(z + h)e~ iuJt }, where the (x,?/)-plane corresponds to the 
undisturbed free surface and z is measured vertically upwards with z = —h the bottom 
of the channel. 

Subsequently ? proved the existence of, and computed the wavenumbers for, the circu- 
lar cross-section case. It should be noted however that experimental evidence for acoustic 
resonances in the case of the circular cylinder is given by ?, pp. 231-232. 
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? provided a theory for determining the trapped-mode frequencies for the thin plate, 
based on a modification of the Wiener-Hopf technique. Further interesting results can 
be found in ? and ?. 

The use of channel Green's functions, developed in § ??, allows the far-field behaviour 
to be computed in an extremely simple manner, whilst the integral equation constructed 
in § ?? enables the trapped modes to be computed in § ??. Conclusions are drawn in § ??. 
Appendix ?? contains the boundary conditions. 



2. Green's functions 

2.1. Construction of equations 

We are concerned with problems for which the solution, 0, is either symmetric or anti- 
symmetric about the centreline of the waveguide, y = 0. The first step is the construction 
of a symmetric and an antisymmetric Green's function, G S (P, Q) and G a (P, Q). Thus we 
require 

(V 2 + k 2 )G s = (V 2 + k 2 )G a = (2.1) 
in the fluid, where V is a gradient operator, 

V-v = 0, V 2 P = V • (v x w). 



In (??) 



G s , G a ~ 1/ (27t) lnr as r = \P - Q\ -» 0, (2.2) 
dG s dG a 



dy dy 
dG s 



on y = d, (2.3) 



= on y = 0, 



8y } (2-4) 

G a = on y = 0, 

and we require G s and G a to behave like outgoing waves as |a;| — > oo. 
One way of constructing G s or G a is to replace (??) and (??) by 

(V 2 + k 2 )G s = (V 2 + k 2 )G a < -5{x - Z)5(y - r?) (2.5) 

and to assume initially that k has a positive imaginary part. 

2.2. Further developments 
Using results from ? we see that this has the integral representation 

/ 7 - 1 [ e - fc Tl^-"l +e- k ^ 2d - y - ri '>}cosk(x-Otdt, 0<y, r, < d, (2.6) 

27t Jo 

where 

' — i(l — * 2 )V2 j t ^ i 
(i 2 -!) 1 / 2 , t>l. 
In order to satisfy (??) we add to this the function 

1 f°° ^coshkj(d-y) , , ^ , 
2nJ w 7sinhfc7d v s; 

which satisfies (??), (??) and obtain 

B{t) = 2e- fe7d coshfc 7 (d-?7). (2.7) 



7(*) = 



